function GAD_Etom_saline_altcontext_crossday_rm_cdf(Data_Summary)
%% WT saline data
G18_saline_RM = [];
for n = 1:size(Data_Summary,2)
    if strcmp(Data_Summary(n).animalName, 'G2_F4_18') && strcmp(Data_Summary(n).drug,'saline')  && isempty(G18_saline_RM)
        [RM, RM_Shuffle] = get_RM_values_ready(Data_Summary(n).Rate_Map_Corr, Data_Summary(n).Shuffled_Rate_Map_Corr_values, Data_Summary(n).RateMap_Corr_Coherent_Rot);
        G18_saline_RM = RM;
        G18_saline_Shuffle = RM_Shuffle;
    elseif strcmp(Data_Summary(n).animalName, 'G2_F4_18') && strcmp(Data_Summary(n).drug,'saline')   && ~isempty(G18_saline_RM)
        [RM, RM_Shuffle] = get_RM_values_ready(Data_Summary(n).Rate_Map_Corr, Data_Summary(n).Shuffled_Rate_Map_Corr_values, Data_Summary(n).RateMap_Corr_Coherent_Rot);
        G18_saline_RM = cat(1,G18_saline_RM,RM);
        G18_saline_Shuffle = cat(1,G18_saline_Shuffle,RM_Shuffle);
    end
end
[G18sf,G18sx,G18slow,G18shigh] = ecdf(G18_saline_RM,'Function','cdf','Alpha',0.05);
[G18ssf,G18ssx,G18sslow,G18sshigh] = ecdf(G18_saline_Shuffle,'Function','cdf','Alpha',0.05);

G28_saline_RM = [];
for n = 1:size(Data_Summary,2)
    if strcmp(Data_Summary(n).animalName, 'G2_F5(2)_28') && strcmp(Data_Summary(n).drug,'saline')  && isempty(G28_saline_RM)
        [RM, RM_Shuffle] = get_RM_values_ready(Data_Summary(n).Rate_Map_Corr, Data_Summary(n).Shuffled_Rate_Map_Corr_values, Data_Summary(n).RateMap_Corr_Coherent_Rot);
        G28_saline_RM = RM;
        G28_saline_Shuffle = RM_Shuffle;
    elseif strcmp(Data_Summary(n).animalName, 'G2_F5(2)_28') && strcmp(Data_Summary(n).drug,'saline')   && ~isempty(G28_saline_RM)
        [RM, RM_Shuffle] = get_RM_values_ready(Data_Summary(n).Rate_Map_Corr, Data_Summary(n).Shuffled_Rate_Map_Corr_values, Data_Summary(n).RateMap_Corr_Coherent_Rot);
        G28_saline_RM = cat(1,G28_saline_RM,RM);
        G28_saline_Shuffle = cat(1,G28_saline_Shuffle,RM_Shuffle);
    end
end
[G28sf,G28sx,G28slow,G28shigh] = ecdf(G28_saline_RM,'Function','cdf','Alpha',0.05);
[G28ssf,G28ssx,G28sslow,G28sshigh] = ecdf(G28_saline_Shuffle,'Function','cdf','Alpha',0.05);

G31_saline_RM = [];
for n = 1:size(Data_Summary,2)
    if strcmp(Data_Summary(n).animalName, 'G2_F5(2)_31') && strcmp(Data_Summary(n).drug,'saline')  && isempty(G31_saline_RM)
        [RM, RM_Shuffle] = get_RM_values_ready(Data_Summary(n).Rate_Map_Corr, Data_Summary(n).Shuffled_Rate_Map_Corr_values, Data_Summary(n).RateMap_Corr_Coherent_Rot);
        G31_saline_RM = RM;
        G31_saline_Shuffle = RM_Shuffle;
    elseif strcmp(Data_Summary(n).animalName, 'G2_F5(2)_31') && strcmp(Data_Summary(n).drug,'saline')   && ~isempty(G31_saline_RM)
        [RM, RM_Shuffle] = get_RM_values_ready(Data_Summary(n).Rate_Map_Corr, Data_Summary(n).Shuffled_Rate_Map_Corr_values, Data_Summary(n).RateMap_Corr_Coherent_Rot);
        G31_saline_RM = cat(1,G31_saline_RM,RM);
        G31_saline_Shuffle = cat(1,G31_saline_Shuffle,RM_Shuffle);
    end
end
[G31sf,G31sx,G31slow,G31shigh] = ecdf(G31_saline_RM,'Function','cdf','Alpha',0.05);
[G31ssf,G31ssx,G31sslow,G31sshigh] = ecdf(G31_saline_Shuffle,'Function','cdf','Alpha',0.05);

G34_saline_RM = [];
for n = 1:size(Data_Summary,2)
    if strcmp(Data_Summary(n).animalName, 'G2_F5(2)_34') && strcmp(Data_Summary(n).drug,'saline')  && isempty(G34_saline_RM)
        [RM, RM_Shuffle] = get_RM_values_ready(Data_Summary(n).Rate_Map_Corr, Data_Summary(n).Shuffled_Rate_Map_Corr_values, Data_Summary(n).RateMap_Corr_Coherent_Rot);
        G34_saline_RM = RM;
        G34_saline_Shuffle = RM_Shuffle;
    elseif strcmp(Data_Summary(n).animalName, 'G2_F5(2)_34') && strcmp(Data_Summary(n).drug,'saline')   && ~isempty(G34_saline_RM)
        [RM, RM_Shuffle] = get_RM_values_ready(Data_Summary(n).Rate_Map_Corr, Data_Summary(n).Shuffled_Rate_Map_Corr_values, Data_Summary(n).RateMap_Corr_Coherent_Rot);
        G34_saline_RM = cat(1,G34_saline_RM,RM);
        G34_saline_Shuffle = cat(1,G34_saline_Shuffle,RM_Shuffle);
    end
end
[G34sf,G34sx,G34slow,G34shigh] = ecdf(G34_saline_RM,'Function','cdf','Alpha',0.05);
[G34ssf,G34ssx,G34sslow,G34sshigh] = ecdf(G34_saline_Shuffle,'Function','cdf','Alpha',0.05);

%% Organize data for WT saline

% get WT saline CDF mean
A = cat(1,G18sf,G28sf,G31sf,G34sf);
WT_saline_F = sort(A);
WT_saline_F = cat(1, 0, WT_saline_F, 1);

% get WT saline CDF x-steps
B = cat(1,G18sx,G28sx,G31sx,G34sx);
WT_saline_X = sort(B);
WT_saline_X = cat(1, -1, WT_saline_X, 1);

% get WT saline CDF high bound mean
C = cat(1,G18shigh,G28shigh,G31shigh,G34shigh);
WT_saline_H = sort(C);

% get WT saline CDF low bound mean
D = cat(1,G18slow,G28slow,G31slow,G34slow);
WT_saline_L = sort(D);

%% MUT saline data
G17_saline_RM = [];
for n = 1:size(Data_Summary,2)
    if strcmp(Data_Summary(n).animalName, 'G2_F4_17') && strcmp(Data_Summary(n).drug,'saline')  && isempty(G17_saline_RM)
        [RM, RM_Shuffle] = get_RM_values_ready(Data_Summary(n).Rate_Map_Corr, Data_Summary(n).Shuffled_Rate_Map_Corr_values, Data_Summary(n).RateMap_Corr_Coherent_Rot);
        G17_saline_RM = RM;
        G17_saline_Shuffle = RM_Shuffle;
    elseif strcmp(Data_Summary(n).animalName, 'G2_F4_17') && strcmp(Data_Summary(n).drug,'saline')   && ~isempty(G17_saline_RM)
        [RM, RM_Shuffle] = get_RM_values_ready(Data_Summary(n).Rate_Map_Corr, Data_Summary(n).Shuffled_Rate_Map_Corr_values, Data_Summary(n).RateMap_Corr_Coherent_Rot);
        G17_saline_RM = cat(1,G17_saline_RM,RM);
        G17_saline_Shuffle = cat(1,G17_saline_Shuffle,RM_Shuffle);
    end
end
[G17sf,G17sx,G17slow,G17shigh] = ecdf(G17_saline_RM,'Function','cdf','Alpha',0.05);
[G17ssf,G17ssx,G17sslow,G17sshigh] = ecdf(G17_saline_Shuffle,'Function','cdf','Alpha',0.05);

G21_saline_RM = [];
for n = 1:size(Data_Summary,2)
    if strcmp(Data_Summary(n).animalName, 'G2_F5(2)_21') && strcmp(Data_Summary(n).drug,'saline')  && isempty(G21_saline_RM)
        [RM, RM_Shuffle] = get_RM_values_ready(Data_Summary(n).Rate_Map_Corr, Data_Summary(n).Shuffled_Rate_Map_Corr_values, Data_Summary(n).RateMap_Corr_Coherent_Rot);
        G21_saline_RM = RM;
        G21_saline_Shuffle = RM_Shuffle;
    elseif strcmp(Data_Summary(n).animalName, 'G2_F5(2)_21') && strcmp(Data_Summary(n).drug,'saline')   && ~isempty(G21_saline_RM)
        [RM, RM_Shuffle] = get_RM_values_ready(Data_Summary(n).Rate_Map_Corr, Data_Summary(n).Shuffled_Rate_Map_Corr_values, Data_Summary(n).RateMap_Corr_Coherent_Rot);
        G21_saline_RM = cat(1,G21_saline_RM,RM);
        G21_saline_Shuffle = cat(1,G21_saline_Shuffle,RM_Shuffle);
    end
end
[G21sf,G21sx,G21slow,G21shigh] = ecdf(G21_saline_RM,'Function','cdf','Alpha',0.05);
[G21ssf,G21ssx,G21sslow,G21sshigh] = ecdf(G21_saline_Shuffle,'Function','cdf','Alpha',0.05);

G45_saline_RM = [];
for n = 1:size(Data_Summary,2)
    if strcmp(Data_Summary(n).animalName, 'G2_F5(2)_45') && strcmp(Data_Summary(n).drug,'saline')  && isempty(G45_saline_RM)
        [RM, RM_Shuffle] = get_RM_values_ready(Data_Summary(n).Rate_Map_Corr, Data_Summary(n).Shuffled_Rate_Map_Corr_values, Data_Summary(n).RateMap_Corr_Coherent_Rot);
        G45_saline_RM = RM;
        G45_saline_Shuffle = RM_Shuffle;
    elseif strcmp(Data_Summary(n).animalName, 'G2_F5(2)_45') && strcmp(Data_Summary(n).drug,'saline')   && ~isempty(G45_saline_RM)
        [RM, RM_Shuffle] = get_RM_values_ready(Data_Summary(n).Rate_Map_Corr, Data_Summary(n).Shuffled_Rate_Map_Corr_values, Data_Summary(n).RateMap_Corr_Coherent_Rot);
        G45_saline_RM = cat(1,G45_saline_RM,RM);
        G45_saline_Shuffle = cat(1,G45_saline_Shuffle,RM_Shuffle);
    end
end
[G45sf,G45sx,G45slow,G45shigh] = ecdf(G45_saline_RM,'Function','cdf','Alpha',0.05);
[G45ssf,G45ssx,G45sslow,G45sshigh] = ecdf(G45_saline_Shuffle,'Function','cdf','Alpha',0.05);

G53_saline_RM = [];
for n = 1:size(Data_Summary,2)
    if strcmp(Data_Summary(n).animalName, 'G2_F5(2)_53') && strcmp(Data_Summary(n).drug,'saline')  && isempty(G53_saline_RM)
        [RM, RM_Shuffle] = get_RM_values_ready(Data_Summary(n).Rate_Map_Corr, Data_Summary(n).Shuffled_Rate_Map_Corr_values, Data_Summary(n).RateMap_Corr_Coherent_Rot);
        G53_saline_RM = RM;
        G53_saline_Shuffle = RM_Shuffle;
    elseif strcmp(Data_Summary(n).animalName, 'G2_F5(2)_53') && strcmp(Data_Summary(n).drug,'saline')   && ~isempty(G53_saline_RM)
        [RM, RM_Shuffle] = get_RM_values_ready(Data_Summary(n).Rate_Map_Corr, Data_Summary(n).Shuffled_Rate_Map_Corr_values, Data_Summary(n).RateMap_Corr_Coherent_Rot);
        G53_saline_RM = cat(1,G53_saline_RM,RM);
        G53_saline_Shuffle = cat(1,G53_saline_Shuffle,RM_Shuffle);
    end
end
[G53sf,G53sx,G53slow,G53shigh] = ecdf(G53_saline_RM,'Function','cdf','Alpha',0.05);
[G53ssf,G53ssx,G53sslow,G53sshigh] = ecdf(G53_saline_Shuffle,'Function','cdf','Alpha',0.05);


%% Organize data for MUT saline

% get MUT saline CDF mean
A = cat(1,G17sf,G21sf,G45sf,G53sf);
MUT_saline_F = sort(A);
MUT_saline_F = cat(1, 0, MUT_saline_F, 1);

% get MUT saline CDF x-steps
B = cat(1,G17sx,G21sx,G45sx,G53sx);
MUT_saline_X = sort(B);
MUT_saline_X = cat(1, -1, MUT_saline_X, 1);

% get MUT saline CDF high bound mean
C = cat(1,G17shigh,G21shigh,G45shigh,G53shigh);
MUT_saline_H = sort(C);

% get MUT saline CDF low bound mean
D = cat(1,G17slow,G21slow,G45slow,G53slow);
MUT_saline_L = sort(D);


%% WT 4h altered contexts data
G28_AB_RM = [];
for n = 1:size(Data_Summary,2)
    if strcmp(Data_Summary(n).animalName, 'G2_F5(2)_28') && strcmp(Data_Summary(n).drug,'AB')  && isempty(G28_AB_RM)
        [RM, RM_Shuffle] = get_RM_values_ready(Data_Summary(n).Rate_Map_Corr, Data_Summary(n).Shuffled_Rate_Map_Corr_values, Data_Summary(n).RateMap_Corr_Coherent_Rot);
        G28_AB_RM = RM;
        G28_AB_Shuffle = RM_Shuffle;
    elseif strcmp(Data_Summary(n).animalName, 'G2_F5(2)_28') && strcmp(Data_Summary(n).drug,'AB')   && ~isempty(G28_AB_RM)
        [RM, RM_Shuffle] = get_RM_values_ready(Data_Summary(n).Rate_Map_Corr, Data_Summary(n).Shuffled_Rate_Map_Corr_values, Data_Summary(n).RateMap_Corr_Coherent_Rot);
        G28_AB_RM = cat(1,G28_AB_RM,RM);
        G28_AB_Shuffle = cat(1,G28_AB_Shuffle,RM_Shuffle);
    end
end
[G28abf,G28abx,G28ablow,G28abhigh] = ecdf(G28_AB_RM,'Function','cdf','Alpha',0.05);
[G28sabf,G28sabx,G28sablow,G28sabhigh] = ecdf(G28_AB_Shuffle,'Function','cdf','Alpha',0.05);

G31_AB_RM = [];
for n = 1:size(Data_Summary,2)
    if strcmp(Data_Summary(n).animalName, 'G2_F5(2)_31') && strcmp(Data_Summary(n).drug,'AB')  && isempty(G31_AB_RM)
        [RM, RM_Shuffle] = get_RM_values_ready(Data_Summary(n).Rate_Map_Corr, Data_Summary(n).Shuffled_Rate_Map_Corr_values, Data_Summary(n).RateMap_Corr_Coherent_Rot);
        G31_AB_RM = RM;
        G31_AB_Shuffle = RM_Shuffle;
    elseif strcmp(Data_Summary(n).animalName, 'G2_F5(2)_31') && strcmp(Data_Summary(n).drug,'AB')   && ~isempty(G31_AB_RM)
        [RM, RM_Shuffle] = get_RM_values_ready(Data_Summary(n).Rate_Map_Corr, Data_Summary(n).Shuffled_Rate_Map_Corr_values, Data_Summary(n).RateMap_Corr_Coherent_Rot);
        G31_AB_RM = cat(1,G31_AB_RM,RM);
        G31_AB_Shuffle = cat(1,G31_AB_Shuffle,RM_Shuffle);
    end
end
[G31abf,G31abx,G31ablow,G31abhigh] = ecdf(G31_AB_RM,'Function','cdf','Alpha',0.05);
[G31sabf,G31sabx,G31sablow,G31sabhigh] = ecdf(G31_AB_Shuffle,'Function','cdf','Alpha',0.05);

G34_AB_RM = [];
for n = 1:size(Data_Summary,2)
    if strcmp(Data_Summary(n).animalName, 'G2_F5(2)_34') && strcmp(Data_Summary(n).drug,'AB')  && isempty(G34_AB_RM)
        [RM, RM_Shuffle] = get_RM_values_ready(Data_Summary(n).Rate_Map_Corr, Data_Summary(n).Shuffled_Rate_Map_Corr_values, Data_Summary(n).RateMap_Corr_Coherent_Rot);
        G34_AB_RM = RM;
        G34_AB_Shuffle = RM_Shuffle;
    elseif strcmp(Data_Summary(n).animalName, 'G2_F5(2)_34') && strcmp(Data_Summary(n).drug,'AB')   && ~isempty(G34_AB_RM)
        [RM, RM_Shuffle] = get_RM_values_ready(Data_Summary(n).Rate_Map_Corr, Data_Summary(n).Shuffled_Rate_Map_Corr_values, Data_Summary(n).RateMap_Corr_Coherent_Rot);
        G34_AB_RM = cat(1,G34_AB_RM,RM);
        G34_AB_Shuffle = cat(1,G34_AB_Shuffle,RM_Shuffle);
    end
end
[G34abf,G34abx,G34ablow,G34abhigh] = ecdf(G34_AB_RM,'Function','cdf','Alpha',0.05);
[G34sabf,G34sabx,G34sablow,G34sabhigh] = ecdf(G34_AB_Shuffle,'Function','cdf','Alpha',0.05);

%% Organize data for WT 4h altered contexts

% get WT saline CDF mean
A = cat(1,G28abf,G31abf,G34abf);
WT_AB_F = sort(A);
WT_AB_F = cat(1, 0, WT_AB_F, 1);

% get WT saline CDF x-steps
B = cat(1,G28abx,G31abx,G34abx);
WT_AB_X = sort(B);
WT_AB_X = cat(1, -1, WT_AB_X, 1);

% get WT saline CDF high bound mean
C = cat(1,G28abhigh,G31abhigh,G34abhigh);
WT_AB_H = sort(C);

% get WT saline CDF low bound mean
D = cat(1,G28ablow,G31ablow,G34ablow);
WT_AB_L = sort(D);

%% MUT 4h altered contexts data
G21_AB_RM = [];
for n = 1:size(Data_Summary,2)
    if strcmp(Data_Summary(n).animalName, 'G2_F5(2)_21') && strcmp(Data_Summary(n).drug,'AB')  && isempty(G21_AB_RM)
        [RM, RM_Shuffle] = get_RM_values_ready(Data_Summary(n).Rate_Map_Corr, Data_Summary(n).Shuffled_Rate_Map_Corr_values, Data_Summary(n).RateMap_Corr_Coherent_Rot);
        G21_AB_RM = RM;
        G21_AB_Shuffle = RM_Shuffle;
    elseif strcmp(Data_Summary(n).animalName, 'G2_F5(2)_21') && strcmp(Data_Summary(n).drug,'AB')   && ~isempty(G21_AB_RM)
        [RM, RM_Shuffle] = get_RM_values_ready(Data_Summary(n).Rate_Map_Corr, Data_Summary(n).Shuffled_Rate_Map_Corr_values, Data_Summary(n).RateMap_Corr_Coherent_Rot);
        G21_AB_RM = cat(1,G21_AB_RM,RM);
        G21_AB_Shuffle = cat(1,G21_AB_Shuffle,RM_Shuffle);
    end
end
[G21abf,G21abx,G21ablow,G21abhigh] = ecdf(G21_AB_RM,'Function','cdf','Alpha',0.05);
[G21sabf,G21sabx,G21sablow,G21sabhigh] = ecdf(G21_AB_Shuffle,'Function','cdf','Alpha',0.05);

G45_AB_RM = [];
for n = 1:size(Data_Summary,2)
    if strcmp(Data_Summary(n).animalName, 'G2_F5(2)_45') && strcmp(Data_Summary(n).drug,'AB')  && isempty(G45_AB_RM)
        [RM, RM_Shuffle] = get_RM_values_ready(Data_Summary(n).Rate_Map_Corr, Data_Summary(n).Shuffled_Rate_Map_Corr_values, Data_Summary(n).RateMap_Corr_Coherent_Rot);
        G45_AB_RM = RM;
        G45_AB_Shuffle = RM_Shuffle;
    elseif strcmp(Data_Summary(n).animalName, 'G2_F5(2)_45') && strcmp(Data_Summary(n).drug,'AB')   && ~isempty(G45_AB_RM)
        [RM, RM_Shuffle] = get_RM_values_ready(Data_Summary(n).Rate_Map_Corr, Data_Summary(n).Shuffled_Rate_Map_Corr_values, Data_Summary(n).RateMap_Corr_Coherent_Rot);
        G45_AB_RM = cat(1,G45_AB_RM,RM);
        G45_AB_Shuffle = cat(1,G45_AB_Shuffle,RM_Shuffle);
    end
end
[G45abf,G45abx,G45ablow,G45abhigh] = ecdf(G45_AB_RM,'Function','cdf','Alpha',0.05);
[G45sabf,G45sabx,G45sablow,G45sabhigh] = ecdf(G45_AB_Shuffle,'Function','cdf','Alpha',0.05);

G53_AB_RM = [];
for n = 1:size(Data_Summary,2)
    if strcmp(Data_Summary(n).animalName, 'G2_F5(2)_53') && strcmp(Data_Summary(n).drug,'AB')  && isempty(G53_AB_RM)
        [RM, RM_Shuffle] = get_RM_values_ready(Data_Summary(n).Rate_Map_Corr, Data_Summary(n).Shuffled_Rate_Map_Corr_values, Data_Summary(n).RateMap_Corr_Coherent_Rot);
        G53_AB_RM = RM;
        G53_AB_Shuffle = RM_Shuffle;
    elseif strcmp(Data_Summary(n).animalName, 'G2_F5(2)_53') && strcmp(Data_Summary(n).drug,'AB')   && ~isempty(G53_AB_RM)
        [RM, RM_Shuffle] = get_RM_values_ready(Data_Summary(n).Rate_Map_Corr, Data_Summary(n).Shuffled_Rate_Map_Corr_values, Data_Summary(n).RateMap_Corr_Coherent_Rot);
        G53_AB_RM = cat(1,G53_AB_RM,RM);
        G53_AB_Shuffle = cat(1,G53_AB_Shuffle,RM_Shuffle);
    end
end
[G53abf,G53abx,G53ablow,G53abhigh] = ecdf(G53_AB_RM,'Function','cdf','Alpha',0.05);
[G53sabf,G53sabx,G53sablow,G53sabhigh] = ecdf(G53_AB_Shuffle,'Function','cdf','Alpha',0.05);


%% Organize data for MUT 4h altered contexts

% get MUT saline CDF mean
A = cat(1,G21abf,G45abf,G53abf);
MUT_AB_F = sort(A);
MUT_AB_F = cat(1, 0, MUT_AB_F, 1);

% get MUT saline CDF x-steps
B = cat(1,G21abx,G45abx,G53abx);
MUT_AB_X = sort(B);
MUT_AB_X = cat(1, -1, MUT_AB_X, 1);

% get MUT saline CDF high bound mean
C = cat(1,G21abhigh,G45abhigh,G53abhigh);
MUT_AB_H = sort(C);

% get MUT saline CDF low bound mean
D = cat(1,G21ablow,G45ablow,G53ablow);
MUT_AB_L = sort(D);

%% Finally, put shuffle data together

% get shuffle CDF mean
A = cat(1,G17ssf,G21ssf,G45ssf,G53ssf,G18ssf,G28ssf,G31ssf,G34ssf);
Shuffle_F = sort(A);
Shuffle_F = cat(1, 0, Shuffle_F, 1);

% get shuffle CDF x-steps
B = cat(1,G17ssx,G21ssx,G45ssx,G53ssx,G18ssx,G28ssx,G31ssx,G34ssx);
Shuffle_X = sort(B);
Shuffle_X = cat(1, -1, Shuffle_X, 1);

% get shuffle CDF high bound mean
C = cat(1,G17sshigh,G21sshigh,G45sshigh,G53sshigh,G18sshigh,G28sshigh,G31sshigh,G34sshigh);
Shuffle_H = sort(C);

% get shuffle CDF low bound mean
D = cat(1,G17sslow,G21sslow,G45sslow,G53sslow,G18sslow,G28sslow,G31sslow,G34sslow);
Shuffle_L = sort(D);


%% plot all data

figure
subplot(1,2,1)
plot(WT_saline_X,WT_saline_F,'Color','#1750AC' ,'LineWidth',3);
hold on
% plot(WT_saline_X,WT_saline_L,'--k','LineWidth',1);
% hold on
% plot(WT_saline_X,WT_saline_H,'--k','LineWidth',1);
% hold on
plot(WT_AB_X,WT_AB_F,'Color','#4DED30','LineWidth',3);
hold on
% plot(WT_AB_X,WT_AB_L,'--k','LineWidth',1);
% hold on
% plot(WT_AB_X,WT_AB_H,'--k','LineWidth',1);
% hold on
% % plot(WT_X_X,WT_X_F,'Color','#4DED30' ,'LineWidth',3);
% % hold on
% % plot(WT_X_X,WT_X_L,'--k','LineWidth',1);
% % hold on
% % plot(WT_X_X,WT_X_H,'--k','LineWidth',1);
% % hold on
plot(Shuffle_X,Shuffle_F,'-k','LineWidth',3);
hold on
% plot(Shuffle_X,Shuffle_L,'--k','LineWidth',1);
% hold on
% plot(Shuffle_X,Shuffle_H,'--k','LineWidth',1);
% hold on
xlabel('RM Correlation Value')
ylabel('Cumulative Fraction of Cells')
xlim([-0.6,1])
set(gca,'FontSize',20)
% set(gca, 'visible', 'off')
hold off

subplot(1,2,2)
plot(MUT_saline_X,MUT_saline_F,'Color','#F53BD6','LineWidth',3);
hold on
% plot(MUT_saline_X,MUT_saline_L,'--k','LineWidth',1);
% hold on
% plot(MUT_saline_X,MUT_saline_H,'--k','LineWidth',1);
% hold on
plot(MUT_AB_X,MUT_AB_F,'Color','#FC6A03' ,'LineWidth',3);
hold on
% plot(MUT_AB_X,MUT_AB_L,'--k','LineWidth',1);
% hold on
% plot(MUT_AB_X,MUT_AB_H,'--k','LineWidth',1);
% hold on
% % plot(MUT_X_X,MUT_X_F,'Color','#4DED30' ,'LineWidth',3);
% % hold on
% % plot(MUT_X_X,MUT_X_L,'--k','LineWidth',1);
% % hold on
% % plot(MUT_X_X,MUT_X_H,'--k','LineWidth',1);
% % hold on
plot(Shuffle_X,Shuffle_F,'-k','LineWidth',3);
hold on
% plot(Shuffle_X,Shuffle_L,'--k','LineWidth',1);
% hold on
% plot(Shuffle_X,Shuffle_H,'--k','LineWidth',1);
% hold on
xlabel('RM Correlation Value')
ylabel('Cumulative Fraction of Cells')
xlim([-0.6,1])
set(gca,'FontSize',20)
% set(gca, 'visible', 'off')
hold off

end